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l.MOF  structured  characterization 
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Figure  SI:  XRD  data  for  HKUST-1  before  and  after  infiltration  with  TCNQ. 

2.  Channel  length  dependence  of  conductivity 

We  measured  the  I-V  characteristics  at  different  temperatures  for  several  devices 
fabricated  with  the  TCNQ@Cu3(BTC)2  thin  film  presented  in  the  main  text.  These  devices 
had  channel  lengths  varying  between  85  and  160  microns.  In  all  cases  we  found  linear  I-V 
characteristics,  and  the  extracted  resistance  as  a  function  of  channel  length  is  shown  in  Fig. 
S2.  The  linear  scaling  with  channel  length  shows  that  the  resistance  is  dominated  by  transport 
in  the  thin  film,  as  opposed  to  contact  resistance. 


Figure  S2:  Resistance  as  a  function  of  the  separation  between  contacts  measured  at  7  different 
temperatures. 

3.  TDTR  measurements  of  thermal  conductivity  and  heat  capacity 

In  TDTR,  the  time  evolution  of  the  surface  temperature  is  measured  through  the 
temperature  dependent  changes  in  reflectivity  (i.e.,  the  thermoreflectance). [1]  In  our  TDTR 
experimental  setup,  laser  pulses  emanate  from  a  Ti:Sapphire  oscillator  with  an  80  MHz 
repetition  rate  and  are  energetically  split  into  pump  and  probe  paths. [21  The  train  of  ultra-short 
pump  pulses  thermally  stimulate  the  A1  metal  dots  (we  confirm  the  thickness  of  each  dot  used 
in  our  TDTR  measurements  with  mechanical  profilometry)  and  time  delayed  probe  pulses 
measure  the  change  in  the  thermoreflectance  of  the  sample  due  to  the  decay  of  the  deposited 
thermal  energy.  Using  lock-in  amplification,  we  analyze  the  ratio  of  the  in-phase  Vm(t)  and 
out-of-phase  Vout(t)  variations  in  the  intensity  of  the  reflected  probe  beam  at  the  modulation 
frequency,/,  of  the  pump  beam  as  a  function  of  delay  time  between  the  pump  and  probe  (up  to 
5.5  ns).  To  ensure  negligible  sensitivity  to  in-plane  transport,  our  pump  and  probe  spot  sizes 
were  focused  to  1/e2  radii  values  of  27.5  and  1 1  pm,  respectively. [1' 3'  41  We  note  that  TDTR 
has  been  used  extensively  to  measure  the  thermal  conductivity  of  thermally  insulative  thin 
films, [5‘9]  validating  our  approach. 


n  21 

TDTR  experimental  details  and  analysis  are  well  documented  in  the  literature,  ’  and 
only  the  pertinent  and  unique  aspects  of  our  specific  experiments  are  discussed  here, 
specifically,  sensitivity  calculations  of  our  measurements  to  the  thermal  conductivity  of  the 
TCNQ@Cu3(BTC)2  MOF  film.  The  A1  metal  dots  were  used  for  TDTR  transduction,  the 
thicknesses  of  which  were  verified  with  profilometry  as  220  nm.  We  analyze  the  TDTR  data 
in  Fig.  4a  of  the  main  document  using  a  model  that  accounts  for  diffusion  through  the  A1  film, 
the  TCNQ@Cu3(BTC)2  MOF  film,  the  100  nm  thermal  oxide  and  the  semi-infinite  Si 
substrate.  To  determine  the  thermal  conductivity  of  the  TCNQ@Cu3(BTC)2  MOF  film,  all 
other  thermophysical  properties  of  each  material  and  interface  in  the  sample  must  be  known. 
Where  heat  capacities  and  thermal  conductivities  of  the  Al,  Si02  and  Si  are  well  documented 
in  the  literature  (including  reduced  thermal  conductivities  of  typical  Al  films  used  in  TDTR 
measurements),  '  "  both  the  heat  capacity  and  thermal  conductivity  of  the 

TCNQ@Cu3(BTC)2  MOF  film  and  thermal  boundary  conductance  across  the  MOF-adjacent 
interfaces  are  unknowns  in  our  analysis.  To  alleviate  this,  we  conduct  TDTR  measurements 
at  a  pump  frequency  of  /  =  1.034  MHz,  which  provides  robust  sensitivity  to  the  thermal 
conductivity  of  the  MOF  (kmof)  while  a  relatively  low  sensitivity  to  the  heat  capacity  of  the 
MOF  (Cmof)  and  thermal  conductances  across  the  Al/MOF  and  MOF/Si02  interfaces 
(^k.ai/mof  and  /?K.Mor/si02-  respectively).  This  is  quantified  in  Figure  S4a,  which  shows  the 
sensitivity  of  the  TDTR  signal  (ratio  of  in-phase  to  out-of-phase  voltages  recorded  by  the 
lock-in  amplifier,  -Vin(t)/-V0ut(t))  to  various  thermophysical  properties  in  the  sample  x,  where  x 
denotes  /cmof,  Cmof,  ^k,ai/mof  or  /tK.MOivsioi-  This  TDTR  sensitivity,  Sx,  is  calculated  by 


Sx(0 


aJ-^1 

L  VoJm 
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where  t  is  the  pump-probe  delay  time.  At  a  pump  frequency  of/=  1.034  MHz,  our  thermal 
signal  is  dominated  by  the  thermal  conductivity  of  the  TCNQ@Cu3(BTC)2  MOF  film.  The 


sensitivity  of  our  TDTR  measurements  at  /  =  1.034  MHz  to  kmof  is  roughly  an  order  of 
magnitude  greater  than  Cmof  over  nearly  the  entire  pump-probe  delay  time. 

Even  with  this  dominant  sensitivity  to  atmof  at/=  1.034  MHz,  we  further  reduce  the 
uncertainty  to  kmof  by  determining  Cmof  through  TDTR  measurements  at/=  8.8  MHz.  By 
changing  the  pump  modulation  frequency,  we  vary  the  depth  under  the  surface  that  is 
thermally  probed  during  a  TDTR  experiment,  and  thereby  vary  the  sensitivity  of  our 
measurement  to  Cmof-[6'  141  As  discussed  above  and  quantified  in  Figure  S3a,  at  1.034  MHz 
our  measurements  are  most  sensitive  to  the  thermal  conductivity  of  the  MOF  more  so  than 
any  other  thermophysical  property  in  the  measurement  volume;  therefore,  at  this  frequency, 
we  fit  our  TDTR  model  to  our  data  by  only  adjusting  the  thermal  conductivity  of  the 
TCNQ@Cu3(BTC)2  MOF,  yielding  relatively  low  uncertainty  in  our  thermal  conductivity 
measurements.  Due  to  the  relatively  low  thermal  conductivity  of  the  MOF,  we  are  insensitive 
to  the  thermal  boundary  conductances  at  each  interface,  and  at  8.8  MHz,  we  are  most  and 
nearly  equally  sensitive  to  both  the  thermal  conductivity  and  heat  capacity  of  the  MOF 
sample,  as  quantified  via  calculations  of  Equation  SI  in  Figure  S4b.  Therefore,  we  use  the 
measured  value  for  thermal  conductivity  at  1.034  MHz  as  input  for  our  analysis  at  8.8  MHz, 
leaving  the  only  free  parameter  for  analysis  at  this  8.8  MHz  as  the  heat  capacity.  Therefore, 
using  these  two  pump  modulation  frequencies,  we  determine  the  thermal  conductivity  and 
heat  capacity  of  the  TCNQ@Cu3(BTC)2  MOF  film  by  iterating  our  TDTR  fitting  analysis 
between  these  two  frequencies.  With  this  dual  frequency  analysis,  we  determine  that  for  our 
TCNQ@Cu3(BTC)2  MOF  film,  Cmof  =  1.20+0.13  MJ  nf3  K1  and  km0 f  =  0.27+  0.04  W  m1 
K  .  The  heat  capacity  obtained  from  our  TDTR  measurements  is  within  25%  of  the  Dulong- 
Petit  (upper)  limit  (1.6  MJ  MJ  m"  K  calculated  using  a  density  of  63.9  mol/F  based  on  1 


TCNQ  per  MOF  unit  cell). 


Figure  S3.  Sensitivities  of  the  TDTR  signal  (ratio  of  in-phase  to  out-of-phase  voltages 
recorded  by  the  lock-in  amplifier,  - V^;n(t)/-  Vout(t))  to  various  thermophysical  properties  in  the 
sample  x,  where  x  denotes  kmof,  CMof,  ^k,ai/mof  or  /?K.M0i7Si02-  These  sensitivities  were 
calculated  via  Eq.  SI  assuming  a  pump  modulation  frequency  of  (a )/=  1.034  MHz  and  (b )/= 
8.80  MHz. 

4.  Molecular  dynamics  simulations  of  thermal  conductivity 

The  Green-Kubo  method  uses  the  auto -correlation  of  equilibrium  heat  flux  J  to 
calculate  the  conductivity  k  from  the  expression 

Kv^m.myt  (2) 

kbi 

where  V  is  the  system  volume,  T  is  the  temperature,  kg  is  the  Boltzmann  constant,  (•)  denotes 

the  canonical  ensemble  average,  and  (J)  =  0  in  equilibrium.  The  microscopic  definition  of  the 
heat  flux  in  terms  of  the  phase  space  of  the  molecular  system  can  be  found  in  Ref.[15]  and 
includes  contributions  from  all  forms  of  interatomic  forces:  covalent  bonds,  short-range 


Lennard-Jones  interactions  and  Coulombic  forces.  The  CVFF1161  and  AMBER1171 
parameterization  were  used  to  represent  the  HKUST  and  TCNQ  molecules,  respectively.  The 
AMBER  parameterization  was  also  used  to  represent  the  intermolecular  forces,  e.g.,  the  Cu-N 
bonding  between  HKUST  and  TCNQ. 

We  employed  2x2x2  super-cells  to  minimize  the  finite  size  effects  based  on  the  work 
of  Huang  et  al.  which  calculated  the  MOF-5  thermal  conductivity1181.  To  ensure  that  the 
canonical  ensemble  was  sampled  sufficiently,  5-10  replicas  derived  from  different  starting 
velocities  and  run  for  1  ns  where  averaged  to  form  the  time -correlation  of  the  heat  flux.  To 
validate  our  calculation  method  and  choice  of  potential,  we  computed  the  thermal 
conductivity  of  MOF-5,  obtaining  k  =  0.25  ±  0.03  W/m-K.  This  is  in  good  agreement  with  the 
value  k  =  0.3 1  ±  0.02  W/m-K  calculated  in  Ref.[18]  and  confirmed  by  experiment  in  Ref [19]. 

We  also  examined  the  dynamics  of  the  system  with  TCNQs  which  interacted  with 
MOF  framework  through  non-bonded  interactions.  Consistent  with  the  enhancement  of  the 
thermal  conductivity  it  appeared  that  the  TCNQs  provided  structured  additions  to  the  MOF 
framework  and  where  also  more  closely  bound  to  one  of  the  two  Cu  pairs  that  comprised  the 
TNCQ  N  to  HKUST  Cu  bond  as  in  Ref.  8. 
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